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CN . Abstract 

The exactly solvable spin-| Ising-Heisenberg model on diamond chain has been considered. We have found the exact results 
for the magnetization by using recursion relation method. The existence of the magnetization plateau has been observed 
at one third of the saturation magnetization in the antiferromagnetic case. Some ground-state properties of the model are 
O ' examined. At low temperatures, the system has two ferrimagnetic (FRIl and FRI2) phases and one paramagnetic (PRM) 
D , phase. Lyapunov exponents for the various values of the exchange parameters and temperatures have been analyzed. It have 
• also been shown that the maximal Lyapunov exponent exhibits plateau. Lyapunov exponents exhibit different behavior for two 
ferromagnetic phases. We have found the existence of the supercritical point for the multi-dimensional rational mapping of the 
spin-i Ising-Heisenberg model on diamond chain for the first time at absence of the external magnetic field and T — > in the 
antiferromagnetic case. 
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I. INTRODUCTION 



The investigation of physical properties of the low-dimensional quantum spin systems with competing interactions 
in an external magnetic field has been a subject of increasingly intense research interest in the recent decades. The 
research interest of these systems has attracted much attention due to the following reasons: first, they can be solved 
exactly by using different mathematical techniques, second, they are realized in the nature, and third, these systems 
present rich thermodynamic behavior, such as the appearance of magnetization plateaus, double peaks structure of 
the specific heat and magnetic susceptibility. 

One of the interesting low-dimensional quantum spin system is the frustrated diamond Heisenberg spin-chain. The 
physical properties of real materials such as copper mineral Cua (003)2(011)2, known as natural azurite (Copper Car- 
bonate Hydroxide) can be well described by using the quantum antiferromagnetic Heisenberg model on a generalized 
diamond chain. The physical properties of the Heisenberg model on diamond chain have been investigated using 
different methods they are full numerical diagonalization and the Lanczos algorithm fl'-'s'l, the decoration-iteration 
transformation 043 j the mapping transformation technique [8], the density- matrix renormalization-group (DMRG) 
and transfer- matrix renormalization-group (TMRG) techniques [9|], Gibbs-Bogoliubov approach [lol |. cluster approach 
the generalized gradient approximations (GGA) [T5|, the density functional theory and state-of-the-art numerical 
many body calculations [T3| . 

Intriguing properties of the azurite made it a good candidate for studying its properties on the low-dimensional 
quantum spin systems. Kikuchi and co-workers |14| have experimentally studied the physical properties of the com- 
pound Cu3(C03)2(OH)2. They have shown that azurite can be regarded as a model substance of a distorted diamond 
chain. The temperature dependence on the magnetic susceptibility and specific heat shows double peak structure 
(around 20 and 5 K) on both magnetic susceptibility and specific heat results. The existence of the magnetization 
plateau at one third of the saturation magnetization has also been experimentally observed in the magnetization 
curve. 

The aim of the recursion relation method is to cut lattice into branches and express the partition function of all 
lattice through the partition function of branches. This procedure will allow to derive one- or multi-dimensional 
mapping for branches of the partition function. After which the thermodynamic quantities of the physical system 
such as magnetization, magnetic susceptibility, specific heat can be expressed through recursion relation. One and 
multi dimensional mapping allows to investigate properties of different models for example Ising model on Husimi 
lattice [15I. [l6l|. zigzag ladde r flTj . triangular lattice two-layer Bethe lattice [l^, [lO], mixed-spin Ising model on a 
decorated Bethe lattice (2ll.l22l|. Q-state Potts model on the Bethe lattice zigzag ladder phase diagrams for 
both ferromagnetic and antiferromagnetic cases, multicritical points, for the spin-1 Ising model on the Bethe lattice 

[Hill- 
In this paper we have investigated some properties of the spin-i Ising- Heisenberg model on diamond chain by using 

dynamical system (recursion relation) approach. Especially we have investigated magnetic properties of the model and 
shown the existence of the magnetization plateau at one third of the saturation value. The investigation of the ground- 
state properties of the model in the A — h plane shows the existence of three phases in the antiferromagnetic case 
and two phases in the ferromagnetic case. Another interesting property of the model has been found by investigating 
the behavior of Lyapunov exponent. Especially we have shown the existence of the plateau in the maximal Lyapunov 
exponent curve. 

The rest of the paper is organized as follows: In the next section using the recursion relation method we derive 
the exact two dimensional recursion relations for the partition function of the spin-i Ising-Heisenberg model on 
diamond chain. The exact results for the magnetization of Ising and Heisenberg spin sublattices have been derived. 
We describe the ground-state properties of the model in A — /i plane. In Sec. IIIII we have discussed the behavior 
of Lyapunov exponent. For the antiferromagnetic case the maximal Lyapunov exponent for the multi-dimensional 
rational mapping is considered and it is shown that near the magnetization plateaus the maximal Lyapunov exponent 
also exhibits plateau structure. The supercritical point at /i = and T ^ has been found. Finally, section IIVI 
contains concluding remarks. 

II. RECURSION RELATION FOR THE ISING-HEISENBERG DIAMOND CHAIN 

Let us consider the spin-i Ising-Heisenberg model on diamond chain with free boundary conditions in the presence 
of an external magnetic field. The Hamiltonian operator of the model is equal to the summation of the plaquette 
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Hamiltonians and can be written as 



N N 
i=l 1=1 

- hH {SI, + SI,) - ^ {^,l + (1) 

where T-Li is Hamiltonian of each plaquette, S"" j, S'^j (a = x, j/, 2:) and /if represent relevant components of Heisenberg 
spin-i and Ising spin-i operators, the parameters J and Ji stand for the interaction between the nearest-neighbouring 
Heisenberg pairs and the nearest-neighbouring Ising and Heisenberg spins, respectively and A is the anisotropy 
parameter. Hamiltonian (1) also includes longitudinal external magnetic fields hn and hi interacting with Heisenberg 
and Ising spins. The first summation in Eq. (1) is corresponding to the interstitial anisotropic Heisenberg spins 
coupling (J and A), the second summation is corresponding to the interaction between the nearest Ising and Heisenberg 
spins and the last two summations are corresponding to the field interaction with Ising and Heisenberg spins. In our 
further calculations we will consider the case when external magnetic field is uniform hn — hi. It is important to 
mention the separable nature of the Ising-type exchange interactions between neighboring Heisenberg dimers which 
are caused from the following commutation rule between different plaquette Hamiltonians: [Hi^T-Lj] = for i ^ j. 
The partition function of the system with Hamiltonian (I) is 

Y exp{-m, (2) 

{tii,Sa,i,Sb,i} 

where /3 = {ksT)^^, kB is Boltzmann constant (hereafter we consider ks — 1) and T is the absolute temperature. 
By cutting diamond chain at Sa,Q and Sb,Q points (central plaquette) into two branches (we denote these branches 
gn{Sa,OT Sb,o) see Fig. 1) the exact recursion relation for the partition function can be derived. After this procedure 
the partition function can be written as 

{Sa.O ,Sb,o} 

where g'^-,{SafiT Sbfi) is contribution of both left and right branches. The sum in Eq. ^ goes over all possible 
combinations of Heisenberg spins Sa.o and Sbfl- Putting into Eq. ([3]) eigenvalues of the operator exp{~ I3[J {S^ o'^'b + 
o^bo + ^^a.o^bo) ~ H^afi + '^'b.o)]} ^au get the partition function expressed through gn{Sa.o, Sb.o) 

Z = e ^+Tg^(n) + e-^+^ff,i( '^^^' ) + e^+^g^( '^^^' )+e ^--ff^(U), (4) 

where by f (up) and J, (down) we denote directions of Heisenberg spins. To find recursion relations for the model we 
need to find relations between gniSafi^ Sbfi) and gn-i{Sa,i, Sb,i)- 

g„{SaQ,SbQ) ^ Y e^'^['^(^».i^^"'.i+'^-i^"^i+'^'^»-i^'-i^ + -''i(^-o+^'-o)^^ (5) 

{Ml:'S„_l,Sb,l} 

* gn-liSa.l, Sb.l)- 

Inserting into Eq. (5) eigenvalues of the operator exp{-P[J{S^ -^S^-^ + S^^^Sl^ + A5^ J + Ji(S'^.o + S^q)pI + 
JifJ-ii^a,! + Sb,i) - HSl^i + 1 -I- 111)]} we can express g„{Safl, Sb,o) through g„-i{Sa,i, Sb,i) 
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AS it can be seen from relations (6) 5n(^^^^) = g„ ( ^^^'^ ) hence our recursion relation will be two-dimensional 
rational mapping. By introducing the following notations 



gn(tt) 

_ gn(U) 
9n[ ^ I 

we can get two-dimensional recursion relation for the partition function 
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Recursion relation (8) plays a crucial role in our further investigation because the thermodynamic quantities like 
magnetization can be expressed through two-dimensional rational mapping. Magnetization for the sublattice of 
Heisenberg spins can be found using the following formula 



= ^ = ^ (9) 

2Z 

In Eq. (9) the sum goes over all possible combinations of 5*0^0 a-nd Sb,o- Putting into Eq. (9) expression for 
the partition function and taking into account the notation (7) we can express magnetization for the sublattice of 
Heisenberg spins through recursion relations which can be written as 

,/A I h n J A h 9 

In the same way we can find magnetization for the sublattice of Ising spins. 

mi =< fli > = < fli >= ■■ . (11) 

In this expression fii is a part of right branch of gniSafi, Sb^)- So to find magnetization for the sublattice of Ising 
spins we need to express gniSa.o, Sb.o) through gn-i{Sa,i, Sb,i)- It is important to mention that this procedure also 
should be done for the partition function. After this procedure the expression for the magnetization of the sublattice 
of Ising spins can be expressed through recursion relation (8): 

■mi =[ifi{xn-i,yn-i) - /2(a;„_i,y„_i))e"^+7'a;„ + (/3(x„„i, - /4(a;„-i, 2/„-i))(e"^+'^ + e^r+w) (12) 

+ ihixn-i,yn^i) ~ /6(a;„-i,2/„-i))e"'i^"^y„]/[(/i(x„_i,y„_i) + /2(a;„_i, y„_i))e"'i^ + ^a;„ 
+ (/3(a;n-i,2/n-i) + /4(a;„-i,2/„_i))(e"^+^T + e^+^^) + (/5(a;n-i, 2/«-i) + /6(a;„-i, y„_i))e"^^"Ty„]^ 
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Expressions (10) and (12) will let us calculate the total single-site magnetization of the spin-i Ising-Heisenberg 
model on diamond chain which can be written as 



mi + 2mii 



(14) 



Figure [2] shows the field behavior of the total magnetization for antiferromagnetic case at the fixed values of interaction 
constants J = 1.5 and Ji = 1, anisotropy parameter A = 1 and different values of the absolute temperature (T). 
At high temperatures the magnetization curve has a monotone structure (Fig. [5] (a)). At lower temperatures the 
plateau of magnetization at one third is arising in magnetization curve (Fig. [2] (b)). Other plots of the magnetization 
curves for the different values of the anisotropy parameter A are displayed in Fig. [2] ((c), (d), (e), (f)). Figures 2 (b), 
(c) and (d) show that the larger positive values of the anisotropy parameter correspond to the larger width of the 
magnetization plateau for the fixed value of the absolute temperature. While for the negative values of anisotropy 
parameter magnetization curves remain the same (Fig. [2](e), (f)). As it can be seen from the figures recursion relation 
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Figure 2: The field dependence of the total magnetization with respect to its saturation value at exchange parameters J = 1.5 
and Ji = 1: (a) T = 0.3, A = 1; (b) T = 0.1, A = 1; (c) T = 0.1, A = 1.5; (d) T = 0.1, A = 0.5; (e) T = 0.1, A = -1; (f) 
T = 0.1, A = -1.5. 



method results are good agrement with other methods results such as the decoration-iteration transformation method 
(for example see [4] figure 3). 

Let us research the ground state of the spin-^ Ising-Heisenberg model on diamond chain via A and h for the 
antiferromagnetic (J = 1.5, Ji = 1) and ferromagnetic (J = —1.5, Ji = —1) models. Depending on the value of ratio 
■j^ and the magnetic field measured in unites of Ji, the system exhibits two ferrimagnetic (FRIl and FRI2) and one 
paramagnetic (PRM) ground-state phases (Fig. 3(a)) for the antiferromagnetic case. Phases FRIl, FRI2 and PRM 
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correspond to the following values of Ising and Heisenberg spins sublattice magnetization: 

FRIl:mi = -0.5,771/^ = 0.5, (15) 
FRI2 : nil = 0.5, 771^ = 0, 
PRM:mi = 0.5,777^ = 0.5. 

Analytically it can be shown that for the fixed values of exchange parameters J = 1.5 and Ji — 1 phase transition 
from FRIl to FRI2 takes place at A = i. Now let us compare the displayed magnetization curves (Fig. [2|) with the 
ground-state phase diagram shown in Fig. 3(a). As it is already mentioned in FRI2 phase the larger positive values 
of the anisotropy parameter (A > i) correspond to the larger width of the magnetization plateau see Fig. 3 (b), (c) 
and (d). In FRIl phase for the fixed values of interaction constants and the absolute temperature the behavior of the 
magnetization curve remains the same (Fig. [2](e), (f)). 

For the ferromagnetic case there are two phases in the phase diagram; the ferrimagnetic (FRI) and paramagnetic 
(PRM) (Fig. 3(b)). Phases FRI and PRM correspond to the following values of Ising and Heisenberg spins sublattice 
magnetization: 

FRI : mi = 0.5, 777^ = 0, (16) 
PRM -.mi = 0.5,777^ = 0.5. 

It can be analytically shown that FRI phase ends on A = — ^ at absence of an external magnetic field. 



III. LYAPUNOV EXPONENT AND SUPERSTABLE POINT 



In this section we will focus on the thermodynamical equilibrium description of the spin-i Ising- Heisenberg model 
on a diamond chain, by studying infinite-size systems. Lyapunov exponents near the magnetization plateau of the 
antiferromagnetic model are interesting to calculate on a diamond chain. It is shown that the behavior of the 
maximal Lyapunov exponent via magnetic field of multi-dimensional rational mapping has a plateau and coincides 
with magnetization one on one-dimensional kagome chain at low temperatures (sl]. It was obtained that the maximal 
Lyapunov exponent had a negative vanishing plateau. 

The following values of Lyapunov exponent can be observed during the investigation. 

1. A < 0. Negative Lyapunov exponents show that the system is dissipative or non-conservative. The systems with 
more negative values of Lyapunov exponent are more stabile. If A = —oo means that we have superstable fixed and 
superstable periodic points. 



7 



2. A = corresponding to neutral fixed point. Zero values of Lyapunov exponents are characteristic for conservative 
systems. At this value of Lyapunov exponents the second-order phase transition takes place. 

3. A > corresponding to unstable and chaotic systems. The systems with positive Lyapunov exponents have 
chaotic behavior. 

In general for the mapping x„ — /(x„_i) Lyapunov exponent X{x) characterizes the exponential divergence of two 
nearby points after n iterations. Lyapunov exponent may be expressed as a limit of mapping stability as [28l - [33 | 

Xix) = lim - In I ^^^^ I . (17) 

n-i-oo n ax 

In multidimensional case with dimension n, exists n of exponents for various directions in space 

n-l 

e^'^jC^^, ...e^" = lim {eigenvalues of the product TT J(a?i))" , (18) 

2 = 

where J(~i) — {tt^) is the Jacobian of the mapping "a^n+i — GCl^n)- For two dimensional mapping (8) we can 
receive the following expression of Lyapunov exponents 

^ n — l 

Ai,A2 — lim —ln{eigenvalues of the product 1 | J{xi,yi)) (19) 

i=0 
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Figure 4: Plot of Lyapunov exponents for spin i Ising-Heisenberg model on diamond chain at exchange parameters (a) the 
antiferromagnetic case at J = 1.5, Ji = 1, A = 1 and temperature T — 0.1 (b) the ferromagnetic case at J = —1.5, Ji = — 1, 
A = 1 and temperature T = 0.1. 



where J{x,y) is the Jacobian of the mapping (8). Expression will let us count up the meanings of Lyapunov 
exponents depending from an external magnetic field (h) at fixed values of constants of interaction (J, Ji), the 
anisotropy parameter (A) and temperature (T). Figure [4] shows the dependence of Lyapunov exponents on an external 
magnetic field for the antiferromagnetic and ferromagnetic cases. As it can be seen from Fig. 2] (a) both values of 
Lyapunov exponent are equal to each other at absence of an external magnetic field. Another interesting property of 
Lyapunov exponent is the existence of the plateau on maximal Lyapunov exponent curve for the antiferromagnetic 
case. It is important to mention that the locations of plateaus on the maximal Lyapunov exponent curve (Fig. |3] (a)) 
coincide with the locations of plateaus on the magnetization curve (Fig. [5] (b)). Such an interesting phenomena of 
Lyapunov exponent has also been observed by investiga ting two, three and six spin exchange interactions Heisenberg 
model on kagome lattice in an external magnetic field [SS] . Figure 2] (b) shows that for the ferromagnetic case the 
maximum value of the maximal Lyapunov exponent tends to zero. 

In Fig. [5] we show the behavior of Lyapunov exponent for the antiferromagnetic case at lower temperature. As it 
can be seen from figure the absolute values of Lyapunov exponents are increasing by decreasing the temperature. 

Next, we turn our attention to the behavior of Lyapunov exponent curves in different ferrimagnetic (FRIl and 
FRI2) phases. For this purpose, Lyapunov exponent curves for fixed values of interaction constants and the absolute 
temperature are plotted in Fig. (6). At low temperatures, the minimum and maximum Lyapunov exponents have 
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different behavior. Only at /i = the maximum value of the minimum Lyapunov exponent equal to value of the 
maximum one, when the system is in FRI2 phase (Fig. 6 (a)), and there is no intersection of Lyapunov exponents for 
h 0. There is a super stable point {Xmax oo), when h = and at T cx) in FRI2 phase. There are two points of 
intersections for the maximum and minimum values of Lyapunov exponents and coincide when h — in FRIl phase 
(Fig. 6 (b)). Lyapunov exponents are tending to zero in thermodynamic limit (T — ^ 0) in FRIl phase at absence of 
an external magnetic field. 

Now let us investigate another interesting property of the recursion relation (8) , namely superstability. First of all 
we will define superstability for the one dimensional recursion relation. Generally one dimensional recursion relation 
= /(a^n-i) is said to be superstable if the following relation takes place |36l - [39| 

where x* is the fixed point of f{x). An other way superstability can be defined by using definition of Lyapunov 
exponent. The system is superstable when 

A = lim i In I ^Cl^ 1= _oo. (21) 

n->oo n dx 

In the same way we can define superstability for two dimensional recursion relations (8). As it is already mentioned 
above for the antiferromagnetic case the absolute values of Lyapunov exponents are increasing by decreasing the 
temperature. Putting values of the exchange parameters ( J = 1.5, Ji = 1 and A = 1) into equation (17) we can see 
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that at thermodynamic hmit at absence of an external magnetic field the following relation takes place for Lyapunov 
exponents 

lim Ai = A2 = -00, (22) 

which shows the existence of the super stable point. 

We have analyzed the behavior of the magnetization for the spin-i Ising-Heisenberg model on diamond chain for 
different values of anisotropy parameter A. For the antiferromagnetic case (J > 0, Ji > 0) at fixed values of exchange 
parameters J = 1.5 and Ji = 1, the temperature T and for A < the magnetization curves have the same appearance 
as in Fig. 2 (e) . The values of characteristic Lyapunov exponent tend to zero at absence of an external magnetic field 
and at T — >■ 0, which means that there is no supercritical behavior for the antiferromagnetic case when A < i. 

For the antiferromagnetic case (J > 0, Ji > 0) at low temperatures and for values of the anisotropy parameter A 
(A > ^) the magnetization curves have the same behavior as shown in Fig. 2 (b). The changes of the anisotropy 
parameter A only brings to the changes of the width of the magnetization plateau at one third. The values of 
characteristic Lyapunov exponent tend to —00 at absence of an external magnetic field and at T — > 0, which means 
that there is a superstable point for the antiferromagnetic case for positive values of anisotropy parameter A. Usually 
a super stable point lies between bifurcation points [S^-IIS- In our case for the spin-i Ising-Heisenberg model on a 
diamond chain there are no bifurcation points but the maximal Lyapunov exponent tends to minus infinity. So we 
get the phase transition in the super stable point at /i = and T — >■ 0. For the first time we get the phase transition 
point at the super stable one. 

IV. CONCLUSION 

By using the recursion relation technique, we have studied magnetic properties of the exactly solvable spin-i 
Ising-Heisenberg model on diamond chain. Recursion relation technique allowed us to construct the exact two- 
dimensional recursion relation for the partition function. The behavior of the total magnetization with respect to 
its saturation value has been investigated. The existence of the magnetization plateau at one third of saturation 
value of magnetization has been observed in the antiferromagnetic case. The ground-state phase diagrams in A — h 
plane show the existence of two ferrimagnetic (FRIl and FRI2) phases and one paramagnetic (PRM) phase in the 
antiferromagnetic case and one ferrimagnetic (FRI) and a paramagnetic (PRM) phases in the ferromagnetic case. 

The properties of Lyapunov exponents were also discussed. The existence of the plateau of the maximal Lyapunov 
exponent curve was observed at low temperatures. It was detected the different behavior for Lyapunov exponent 
curves in two ferrimagnetic phases. We have shown that for the antiferromagnetic case in the thermodynamic limit 
(T — >■ 0) both values of Lyapunov exponent tend to — cx3 at absence of an external magnetic field which is corresponding 
to the superstable point. 
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